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Abstract. In this note we present a Hopf algebra description of a bosonic quantum 
^ model, using the elementary combinatorial elements of Bell and Stirling numbers. 

Our objective in doing this is the following. Recent studies have revealed that 
perturbative Quantum Field Theory (pQFT) displays an astonishing interplay between 
analysis (Ricmann Zeta functions), topology (Knot theory), combinatorial graph 
j«^ theory (Feynman Diagrams) and algebra (Hopf structure) . Since pQFT is an inherently 

f^ complicated study, thus far not exactly solvable and replete with divergences, the 

^N essential simplicity of the relationships between these areas can be somewhat obscured. 

. . The intention here is to display some of pevious-mentioned structures in the context of 

^ ^ a simple bosonic quantum theory; i.e. a quantum theory of non-commuting operators 

S^ which do not depend on spacetime. The combinatorial properties of these boson 

S-H creation and annihilation operators, which is our chosen example, may be described by 

graphs, analogous to the Feynman diagrams of pQFT, which we show possess a Hopf 
algebra structure. Our approach is based on the quantum canonical partition function 
for a boson gas. 
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1. Bell and Stirling numbers 

A basic combinatorial construct is the Bell number, B{n). This gives the number of 
ways in which n distinct objects may be distributed among n identical containers, some 
of which may remain empty. The first few values of B (n) are 

B{n) = 1, 1, 2, 5, 15, 52, 203, . . . , n = 0, 1, 2, . . . . 

These values grow rapidly, but less so than n! 

Related to the Bell numbers are the Stirling numbers (of the second kind) S{n, k), 
which are defined as the number of ways of putting n different objects into k identical 
containers, leaving none empty. From the definition we have B{n) = I]fc=i S{n, k). 

A compact way of defining a combinatorial sequence is by means of a generating 
function. For the Bell numbers, the exponential generating function is given by 

Ei?H- = exp(e--l) (1) 

n=0 ""• 

A slight generalization of the exponential generating function for the Bell numbers 
is given by that defining the Bell polynomials Bn{y): 

oo n 

exp{y{e^-l)) = j:Bniy)^ (2) 

Note that 



n=o ^■ 



Bn{y) = Y.S{n,k)y'. (3) 

k=0 

1.1. Normal order 

Although somewhat unfamiliar to physicists, the Bell and Stirling numbers are 
fundamental in quantum theory. This is because they arise naturally in the normal 
ordering problem of bosonic operators. By the normal ordering operation M f{a, a'^) we 
mean reorder the boson operators in /(a, a^) so that all annihilation operators are on 
the right. For canonical boson creation and annihilation operators {a\ a] satisfying 



a, a 



1, the Stirling numbers of the second kind S{n, k) intervene through [U |2] 

U{a^aY = Y,S{n,k){a^faK (4) 



fc=i 



The corresponding Bell numbers B{n) = J2k=i ^i^y ^) ^^^ simply the expectation values 

B{n) = {z\{a^ay^\z),=, (5) 

taken in the coherent state defined by 

a\z) = z\z) (6) 

for z = 1. In fact, for physicists, these equations may be taken as the definitions of the 
Stirling and Bell numbers. 
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1.2. Graphs 

We now give a graphical representation of the Bell numbers, based on work of Brody, 
Bender and Meister [3l H] , which we have extended in [5] . Consider labelled lines which 
emanate from a white dot, the origin, and finish on a black dot, the vertex. We shall 
allow only one line from each white dot but impose no limit on the number of lines ending 
on a black dot. Clearly this simulates the definition of S{n, k) and B{n), with the white 
dots playing the role of the distinguishable objects, whence the lines are labelled, and 
the black dots that of the indistinguishable containers. The identification of the graphs 
for 1, 2 and 3 lines is given in Figure 1 below[|} We have concentrated on the Bell number 

n=1: O • V, 

S(1,1)=1 B(1)=1 



O • 

n=2: V,' 

O • 

S(2,2)=1 S(2,1)=1 B(2)=2 




O 




n=3: O • v/ 

O • O • V, 

S(3,3)=1 S(3,2)=3 S(3,1)=1 B(3)=5 

Figure 1. Graphs for B(n), n = 1, 2, 3. 

sequence and its associated graphs since, as we shall show, there is a sense in which this 
sequence of graphs is generic in the evaluation of the Quantum Partition Function. By 
this we mean the following: We first show below that the Bell number sequence arises 
from the Partition Function of a non-interacting boson model. Then, when interactions 
are present these may be incorporated by the use of suitable strengths associated with 
the vertices of the graphs. That is, we can represent the combinatorial sequence of an 
interacting model by the same sequence of graphs as in Figure 1, with suitable vertex 
multipliers (denoted by the V terms in the same figure). 

2. Partition Function Integrand 

We now show the relation of the preceding considerations to the computation of 
the Canonical partition Function Z in Quantum Statistical Mechanics, and introduce 

I The diagrams of our Figure 1 correspond to those of Figure 3 of reference [3] , with an interchange of 
black spots and white spots. We have regrouped them to make the relation with B(n) and S(n.k) more 
transparent. 
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the Partition Function Integrand, which is directly related to the combinatorial Bell 
numbers. 

2.1. Free boson gas and Bell polynomials 

The canonical partition function associated with the hamiltonian H is given by 

Z = TTexp{-/3H). (7) 

We take the elementary case of the hamiltonian for the single-mode free boson gas 
H = ea'^a (ignoring an additive constant), e > 0. The usual computation of the partition 
function, exploiting the completeness property I]^o l^)(^l — -^' i^ immediate: 

Z = Trexp(— /3ea^a) (8) 

oo 

= Y. {n\e-^'''\n) (9) 

n=0 

oo 

= 5]e-^^" (10) 

n=0 

l-e-^'Y\ (11) 



However, we may use any complete set to perform the trace. We choose coherent states 
as defined in Eq.([6]) above, which are explicitly given by 

|;2) = e-l^l'/2^(z"/n!)at"|0). (12) 

n 

For these states the completeness or resolution of unity property is 

- fdh\z){z\=I= fdfi{z)\z){z\. (13) 

The appropriate trace calculation is 

Z = - f d^z{z\exp{-(3ea^a)\z) = (14) 

TT J 

= - fd^z{z\ : exp(a^a(e-^' - 1)) : \z) , (15) 

n J 

where we have used the following well-known relation [6l [7] for the forgetful normal 

ordering operator : /(a, a^) : which means "normally order the creation and annihilation 

operators in / forgetting the commutation relation [a, a^] = 1"r| 

M exp{xa^ a) =: exp(a^a(e^ — 1)) : . (16) 

We therefore obtain, integrating over the angle variable 9 and the radial variable r = \z\, 

T /'27r /-oo , , 

Z = - dO rcirexp(r2(e-''^-l)) , (17) 

TT JO JO ^ ^ 

which gives us Z = (1 — e^^^)~ as before. 

§ Of course, this procedure may alter the value of the operator to which it is applied. 
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We rewrite the above equation to show the connection with our previously-defined 



combinatorial numbers. Writing y = r and x = — /3e, Eq.(17) becomes 



Z = |rfyexp(y(e^-1)). (18) 



This is an integral over the classical exponential generating function for the Bell 
polynomials as given in Eq.(|2]). This leads to the combinatorial form for the partition 
function 



X 



dyJ2Bn{y)^. (19) 



n=0 ^■ 



Although Eq.(19) is remarkably simple in form, it is often by no means a 
straightforward matter to evaluate the analogous integral for other than the free 
boson system considered here. Further, it is also clear that we may not interchange 
the integral and the summation, as each individual y integral diverges. We shall 
therefore concentrate in what follows on the partition function integrand (PFI) F{z) = 
{z\exp{—f3H)\z), whence Z = J F{z) dfi{z), to give a graphical description of a 
perturbation approach. The function F maps coherent states \z) to (real) numbers, 
and is therefore a functional on the coherent states.. 

2.2. General partition functions 

We now apply this graphical approach to the general partition function in second 
quantized form. With the usual definition for the partition function Eq.([7]). In general 
the hamiltonian is given hj H = ew{a, a^), where e is the energy scale, and w is a string 
(= sum of products of positive powers) of boson creation and annihilation operators. 
The partition function integrand F for which we seek to give a graphical expansion, is 



Z{x) = / F{x,z)di^{z), (20) 

where 

F{x,z) = {z\ exp{xw)\z) = {x = —f3e) 

oo n 

= j:{^\w-\z)^ 

= E w.{z) ^ 

n=0 ""• 

/ oo n\ 

= exp EK(^)-r , (21) 



n=l 



n\ 



with obvious definitions of Wn and V^. The sequences {Wn} and {Vn} may each be 
recursively obtained from the other [8]. This relates the sequence of multipliers {Vn} 
of Figure 1 to the hamiltonian of Eq.([7]). The lower limit 1 in the Vn summation is a 
consequence of the normalization of the coherent state | z) . 

At this point the analogue of this model to perturbative quantum field theory 
becomes more transparent. For example, to calculate the partition Function Z{x), we 
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must integrate over the Partition Function Integrand F{x, z). Here, the integration over 
the coherent state parameter z plays the role of the space-time integration of Quantum 



Field Theory. As noted above with reference to Eq.(19) for the non-interacting case, 
and as in pQFT, term-by-term integration, corresponding to fallaciously interchanging 
the summation and integration actions, results in infinities at each term, as can be 
directly verified even in the non-interacting case. Thus even here an elementary form of 
renormalization is necessary, if we have to resort to term-wise integration. 

Considerations such as these have led many authors, as noted in reference [T0|, to 
consider a global, algebraic approach to pQFT; and this we now do in the context of 
our simple model. 

3. Hopf Algebra structure 

To describe the Hopf Algebra structure of our model, which we shall refer to as BELL 
below, we first introduce a basic Hopf Algebra, which we call POLY, generated by a 
single parameter x. This is essentially a standard one- variable polynomial algebra on 
which we impose the Hopf operation of co-product, co-unit and antipode, as a useful 
pedagogical device for describing their properties. The Hopf structure associated with 
our model described above is simply a multi- variable extension of POLY. 

3.1. The Hopf Algebra POLY 

POLY consists of polynomials in x (say, over the real field IZ, for example). The standard 
algebra structure of addition and associative multiphcation is obtained in the usual way, 
by polynomial addition and multiplication. The additional Hopf operations are: 

(i) The coproduct A : POLY — > POLY x POLY is defined by 

A(e) = e X e (unit e) 

A(a;) =xxe-|-exa; (generator x) 

A{AB) = A{A)A{B) otherwise 

so that A is an algebra homomorphism. 

(ii) The co-unit e satisfies e(e) = 1 otherwise e{A) = 0. 

(iii) The antipode S : POLY — )• POLY satisfies S{e) = e; on the generator x, 
S{x) = —X. It is an anti-homomorphism, i.e. S{AB) = S{B)S{A). 

It may be shown that the foregoing structure POLY satisfies the axioms of a 
commutative, co-commutative Hopf algebra. 

3.2. BELL 

We now briefly describe the Hopf algebra BELL which is appropriate for the diagram 
structure introduced in this note, defined by the diagrams of Figure 1. 
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(i) Each distinct diagram is an individual basis element of BELL; thus the dimension 
is infinite. (Visualise each diagram in a "box".) The sum of two diagrams is simply 
the two boxes containing the diagrams. Scalar multiples are formal; for example, 
they may be provided by the V coefficients. 

(ii) The identity element e is the empty diagram (an empty box). 

(iii) Multiplication is the juxtaposition of two diagrams within the same "box" . BELL 
is generated by the connected diagrams; this is a consequence of the Connected 
Graph Theorem [S]. Since we have not here specified an order for the juxtaposition, 
multiplication is commutative. 

(iv) The coproduct A : BELL — > BELL x BELL is defined by 

A(e) = e X e (unit e) 
^(y) =yy<e + exy (generator y) 
A{AB) = A{A)A{B) otherwise 
so that A is an algebra homomorphism. 
(v) The co-unit e satisfies e(e) = 1 otherwise e{A) = 0. 

(vi) The antipode S : BELL — > BELL satisfies S{e) = e; on a generator x, 
S{y) = —y. It is an anti-homomorphism, i.e. S{AB) = S{B)S{A). 

3.3. BELL as an extension o/POLY 

It may be seen that BELL is a multivariable version of POLY. To show this, we code 
the diagrams by letters. We use an infinite alphabet Y = {yi,y2,y3 ■ ■ ■} = {yk}k>i and 
code each connected diagram with one black spot and k white spots with the letter y^. 
An unconnected diagram will be coded by the product of its letters. 
In this way, the diagrams of Figure 1 are coded as follows: 

• first line : yi 

• second line : yf and ?/2 

• third line : yf and yiy2 and y^ . 

Thus one sees that each diagram of weight n with k connected components is coded 
bijectively by a monomial of weight n (the weight of a monomial yi^yi^ ■ ■ ■ Z/v is just the 
sum of the indices Z]j=i ^j) and k letters. The algebra BELL is coded by commutative 
polynomials in the infinite alphabet Y; that is, the coding is an isomorphism BELL — )■ 
7^[y]. As an aside, one may note that the basis elements of BELL are sometimes 
referred to as forests. As above, it may be shown that the foregoing structure BELL 
satisfies the axioms of a commutative, co-commutative Hopf algebra. 

4. Discussion 

The objective of this note was to introduce a very simple Hopf structure associated with 
standard bosonic creation and annihilation operators, in the context of the evaluation 
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of the canonical partition function of quantum statistical mechanics. We did this via a 
diagrammatic description applied to a non-interacting boson gas, but implied that the 
algebraic description was general enough to include interactions as scalar coefficients. 
The inspiration for this task arose from considering recent work on perturbative quantum 
field theory (pQFT) [TU] , where it has been shown that a Hopf description is available 
for this far more complicated system. It is instructive to consider a much more 
straightforward system, such as the one treated here, where the operators do not depend 
on space or time. Nevertheless, we have shown that even such a basic system does exhibit 
some of the features of the more comphcated case, in particular the structure of a Hopf 
algebra, which we called BELL. This may be thought of as a simple solvable model in 
its own right. 

However, one may also ask wherein does this simple structure sit within the full 
pQFT structure? The strategy which we adopt, and which we describe in further work, 
is to generalize the algebraic structure, and thereby produce Hopf algebras of sufficient 
complexity to emulate those associated with pQFT. In subsequent work we show that 
by suitable deformation of the Hopf algebra described here we obtain structures related 
to those arising from more realistic models, such as those associated with pQFT. 
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